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ABSTRACT 


We present a theoretical model for the plasma in the Jovian 
magnetosphere whose pressure is comparable to the corotational 
energy density. The model predicts a thin current sheet of 
1 R, - 2 R_ half-thickness , The current sheet lies almost pre- 

U U 

<5isely in the magnetic equatorial plane and is not appreciably 
warped as has been siaggested previously. 


I. INTRODUCTION 


Observations of the magnetic field structure and energetic 
particle fluxes in the Jovian magnetosphere have been interpreted 
by many authors as indicative of thin disk-like magnetosphere [e.g. , 

Smith et al. , 197U{ Van Allen et al. , 197^> Goertz et al. , 19753* 

This topology was first suggested by Gledhill [19673* The Jovian 

magnetosphere contains a thin current sheet which is presumably due 

to the rapid rotation of the planet. It has been proposed by Van 

Allen et al. [1974] and by Goertz et al. [19753 that the current 

sheet lies almost precisely parallel to the magnetic equatorial plane 

which is inclined with respect to the rotational equatorial plane by 

an angle o: w 10°. It has been suggested that this may not be the 

case for Jupiter because the centrifugal force will warp the current 

sheet so that it becomes parallel to the rotational equator [see, 

e.g.. Hill et al. , 1974]. We will show below that this is not true 

for the case of Jupiter where the plasma pressure in the current 

sheet is comparable to or larger than the corotational energy density. * 

(This possibility was also mentioned by Hill et al. [1974] but not 

studied in detail*) In fact, the deviation of the current sheet 

from the magnetic equatorial plane may be very small (< 1 Rj) even 

at a distance of 100 Rj, 

In section II we will review some observations relevant to 
the problem. In section III we will develop a model illustrating 
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the effects of rotation on the Jovian magnetosphere • This model Is 
very similar to that of Axford and Gleeson [1975 1 "but differs in so 
far as it allows for the temperature of the plasma to be large whereas 
Axford and Gleeson restrict their treatment to the case of loW"energy 
plasma. 


II. OBSEHVATIONS 


By now it seems to be an accepted conclusion that the Jovian 
magnetosphere corotates with the planet. Due to this corotation 
and thus wobbling of the magnetosphere (due to the tilt of the dipole) 
the two spacecrafts Pioneer 10 and 11 sampled different regions of 
latitude as the planet rotated beneath the spacecraft. The ener- 
getic particles seem to be confined to the vicinity of a current 
sheet (or neutral sheet). Also, the magnetic field is strongly 
distended which is indicative of a strong azimuthal current confined 
to a narrow region in latitude. The current sheet is presumably 
neither precisely parallel to the magnetic (dipole) equator nor 
precisely parallel to the rotational equator. It lies somewhere 
between these two planes. Figures la and lb show two models for 
the shape of the current sheet. Figure la depicts a warped current 
sheet tending to become parallel to the rotational plane [Smith et al . , 
197 ^; Hill et al. , 1974]. Figure lb shows an essentially rigid sheet 
aligned with the magnetic equatorial plane [Van Allen et a" - . '741. 

There are two observations which may help to decide which configura- 
tion is actually realized in the Jovian magnetosphere. 

1. Smith et al. [1974] calculate the angle between the 
perturbation field b (observed field minus dipole field) and the 
axis of rotation 1k as well as the angle between b and the magnetic 
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dipole axis iS# They find that eos”^ (1^ • fi) does not vary appre- 
ciably along the trajectoxy of Pioneer 10 except for short times when 
the spacecraft presumab.iy dips into the current sheet and encounters 
a rapidly varying field. The angle cos"^ (li • however, shows a 
more or less sinusoidal variation except for the short times when 
the spacecraft dips into the current sheet. The angle cos (R • b) 
tends to approach 90 “ as the spacecraft moves away from the planet. 
Smith et al. conclude correctly that the perturbation field in the 
limited region covered by Pioneer 10 becomes parallel to the rotational 
equator. But they further conclude, not necessarily correctly, that 
this proves that the current sheet itself becomes parallel to the 
rotational equator. 

We think that this conclusion is not necessarily correct for 
the following reason: As the current sheet is rotated around the 
planet it also moves up and down (see Figure la). Thus the perpen- 
dicular distance z of the spacecraft frwn the current sheet varies 
in time. A constant angle cos**^ would hence ln^jly that the 

dliection of b does not change with perpendicular distance z. In 
other words 




/b \ b 

1 3b 

3b 

z 1 r 

~r 7 

b — . 
1 r az 

' ^ 

z az 


= 0 




where b = b z + b r and z is either along R or M. 
z r ^ 
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At large distances iiuiii the current sheet the current density 
5b^/Sz is essentially zero and the angle cos"^ (R » fi) should only 
be a constant if 





0 . 


Thus the fact that the cos”^ (R . b) is a constant is consistent 
with the idea of a thin current sheet parallel to the rotational axis 
only if bp decreases as l/r. This, however, does not seem to be the 
case (see the magnetic field data of Smith et al. [197^]). 

On the contrary, it seems much more plausible that b^/b 
should vary with z in some regular way. Thus the fact that 
cos ^ (R . ^) is a constant and cos”^ (l3 . b) is not, does neither 
prove that the current sheet is parallel to the rotational axis 
nor parallel to the magnetic axis. We conclude that the observation 
of Smith et al. is inconclusive as to whether or not the current 
sheet is warped. This con<’lusion is supported by the fact that the 
model field of Goertz et al,, which is based on a current sheet 
precisely parallel to the magnetic equatorial plane, reproduces the 
observations of Smith et al. quite satisfactorily. Figure 2 dis- 
plays the angles cos" (R . b) and cos" (M . b) calculated from the 
model; the dashed lines shows the observations of Smith et al. 

2. The second observation is already contained in the first. 
From the drops in magnetic field strength and field direction. 
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Smith et al. concluded that Pioneer 10 dipped into the current ahoot 
on many occasions (e.g., when it was at 86 Rj during the outbound 
pass) actually passed through the current sheet, A passage through 
the current sheet is only possible if the Joviographic latitude of 
the current sheet is larger than the spacecraft's latitude. At 
86 Rj the latitude of PIO was 8.6®, Thus the current sheet must 
have been at a latitude of at least 8.6®. Or, the magnetic lati- 
tude of the current sheet must have been leas than a - 8,6® ~ 1.4®. 
As Figiure la clearly shows, any appreciable warping of the current 
sheet will increase its magnetic latitude above this value. From 
the fact that Pioneer 10 actually crossed the current sheet several 
times during the outbound pass we conclude that the warping must 
everywhere be less than about 1®, i.e., nowhere near the 8® implied 
by the warped model of Figure la. 

Finally we point out thau the model field of Goertz et al. 
reproduces not only the observed directions of ^ but also the mag- 
nitude of the observed fields to better than 5^ accuracy (figure 5). 
(it should, however, be noted that their model applies only to the 
outbound pass of Pioneer 10, i.e., to Jupiter's magnetotail.) Al- 
though the good agreement between the model field and the observa- 
tions is not a proof that the current sheet lies in the magnetic 
equatorial plane, it tends to support that assumption. 

In conclusion we would like to say that no observation 
clearly distinguishes between the two models, although the rigid 
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disk nodtil of Van Aller seems more in agreement with the data than 
the warped model. In the following section we will show that theory 
predicts very little warping. 
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III. OHE MODEL 

Consider a cylindrical coordinate system (r, cp, z) with the 
z axis parallel to the magnetic dipole axis. Let '^he whole magneto- 
sphere rotate with an angular velocity n about an axis inclined at an 
angle a with respect to the z axis. We assume <x small but not zero. 
Initially we assume that a thin plasma sheet lies in the plane z = 0. 
It is thin enough so that its magnetic field is equivalent to that 
of an azimuthsl current sheet. We also assume that only meridional 
field components exist, i. e. , B « 0. The last two assumptions 
are identical to those made by Gleeson and Axford. Following them 
we also assume that close to the current sheet 


i-« 1- 

dr dz • 


Thus V ' B = 0 reduces to 



B « B (r) 
z zo' ' 


( 1 ) 


(n) 


Then: 
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The fort e equations for the plasma sheet in equilibrium are then 


n 1 

-pf2‘ (r COG I + z nin 0 cos + — JB - — r?- cos 

f)r c cp z >’ 
r 

(4) 


j p. , QMtP 

-p!2' (r coo + z sin t) oin oi = - — ^ ~ j ” 


B - — rr- sin 


dz c <p r ^r’ 


tan i'’- = z/r , 


(5) 

(6) 


where p is the placma density, p is pressure and GM_ = 1.:’6 y x 10 cm /s 

O ' 

[Anderson et al. , 1974]. We now restrict ourselves to small values of f'., 

l.e., regions close to the magnetic equator. We also adopt 

sin « cos (X, Then near the equator (z « r) the equations reduce to 


-pH' r 


^ + 1 j 

9r c '^(p 


CM. 

U 

— P 
r 


(7) 


-pil' r cos (i sin B 

3z c cp r 


( 8 ) 




We know that 


B=B° + h B=B° + b 
r r r z z z 


(9) 


L' 


where t! e euperBcript D refera to the dipole field. We now write; 


p k(r,z) . (10) 

We recover the result of Gleeson and Axi'ord if k « 1. However, if 
the magnet ocphe re plasma is supplied by the Jovian ionosphere (as in 
the models of loannldis and Brice [l97l] and Goertz [1975 1) k > 1. 

Also Frank et al. ll?i'5] find thermal plasma in the Jovian magneto- 

P P 

sphere with p *=» 10 . Thus k « 1 does not seem to be a valid 

assumption for the Jovian magnetosphere. 

Equation (8) can be rewritten as 


-pli r cos 


(z sin a + 




O 



( 11 ) 


In the symmetry plane of the current sheet the right-hand ride of 
this equation is identically equal to zero, because both tue kijictic 
pressure p and magnetic pressure b‘^8r' are symmetric about the current 
sheet. Since the first term of the left-hand side is non-zero we 
must have a deviation of the current sheet from the plane z = 0, 
where = 0. Let us assume that the current sheet is symmetric 
about the surface z = d(r). We require then that 
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1^- 


? 1 9^ r 

p ^ r cos a sin a = S. ' 

J+ir az rlz=Td(r) 


( 12 ) 


Foi Jupiter 


12zr r_i 


(«hera z anc r are expreaaed ia unite of Jovian radii), m the 

current eheet ue then have, aeeuming the left-hand aide of Eq. (u) 
to be zero^ 


p + = c(r) . 


( 15 ) 


G«rte et al. [1975] have ahc™ that eloae to the ma«netio equator 
( z « r) aud for distances r > . o 




a = 0.7 , 


P « 1 Rj . 


( 14 ) 


We also require that p -4 o for 


z -♦ ». Thus 


p = 


1 o 2 

UJ 


( 15 ) 
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Insert in,; Eqs, (l4) and (15) into Eq. (12) yields 
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D cos JK sin a b_ 
z 


D "5 

where B r"^ » iKJ. 
z 

Combining Eqs. (7), (lO), (l4), and ( 15 ) gives 
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zo 


k(a + 1 ) 4 (1 


GMj 
n r 


\1 


(16) 


(17) 


Goertz et al. (1975) have shown that the perturbation field close 
to the magnetic equator (z « r) has the form 


acb D 
0 

a+^ 


( 18 ) 


where c =- 10. Since B = B^ 4 - b we find 

zo z z 


I, = 


(19) 


^o ack 4- k(a 4 - 1) 4 - 


1 - ' 

fi r i 


With the values of Goertz et al. (1975) for ihe constants a, c, and 

b we find 
o 


Vn 
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1 

k 


II - 


GM_ 



( 20 ) 


2 

where the inequality holds for r > S, i.e. for 1> GMj/fi r". Equa- 
tion (2o) tells us that even for a high-energy plasma (k > 1) the 
current sheet is thin. In fact, for k » «, D VEu:‘ies from D r- 1.6 
at r = 20 to D ss 1 at r a lOO. This is in excellent agreement with 
the value of D = 1 that was assumed in the magnetic field model of 
Goertz et al. at r ■> 20. 

Finally we calculate the deviation of the current sheet from 
the magnetic equatorial plane by inserting equation (19) into 
equation (l6). For d « r, which should be the case for r > 10, 
we find 


d « 1 cos o. sin a 

? 5 


k(l + a + ac) + jl - 


GMj 

0 r 


( 21 ) 


We recover the result of Hill et al. (1974) for k = 0 and 
1 > GMj/n^r^. 1 > GMj/n^r^ was also asr'jmed by Hill et al. (19T4). 
It is certainly true in the regions where equation (21) is valid, 
namely r s 20. For k » 1 the deviation of the current sheet from 
the magnetic equatorial plane is about 10 times smaller than the 
value predicted by Hill et al. (1974), It should be noted that 
k > 1 corresponds precisely to the case, mentioned by Hill et al., 


where th ! magnetic mirror force equals or exceeds the* centrifugal 
stresses . 


We are now in the position to calculate the density in the 
current sheet. By using equations (lo) and (15) we find 


Po 




b 

0 


, 2 
4Trn r 




( 22 ) 


The plasma density for a hydrogen plasma is then 

N w 3 X 10^ 1 

° k 5.'V 

r 


(23) 


The results of Frank et al. (1975) indicate a value for k ~ 5 at 
r a- 10 which is, however, uncertain to at least a factor of 2. 
Equation (23) predicts a value of N 10 /k cm'*'"^ at r = 10, which 
for k =- 5 is an order of magnitude larger than the plasma density 
observed by Frank et al. We regard this discrepancy not as too 
serious for two reasons; i) equation (23) is valid only for r > 20, 
where b and b are well described by equations ( 19 ) and (18), 
ii) the temperatures and densities reported by Frank et al. may well 
be tincertain by factors of 2-5 each. 

Finally we note that the self-consistent current density 
calculated from the drift velocities 
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■4 




-* 4 

F X B 


is numerically conoictent with the value of 


j=_c_vxB = _c_ 

4n 4 tt 

to within a factor of 1.5. The slight discrepancy is purely due to 
uncertainties in the numerical values of the parameters a and c and 
the neglect of the current density due to ab^/ar. 


IV. DISCUSSION 


The treatment of the current sheet above has two inailor 
shortcomlnRo , 

(i) It la baaed on a magnetic field model which describes 
well only the mnffnetic field obaervations obtained on the outbound 
lefi of Pioneer 10. It is by no means clear that a similar analysis 
can be performed for the front side of Jupiter's magnetosphere. 
Clearly, on analysis of the magnetic field data obtained by 
Pioneer 10 and 11 while they traversed the front side magnetosphere, 
similar to that by Ooertz et al. (1975) is needed, 

(ii) . The solutions had to be left in terms of on unspeci- 
fied parameter k. Unfortunately no plasma measurements are avail- 
able in those regions where the analysis is valid, k must thus be 
determined from some kind of model calculation. Quantitative models 
for the plasma terapei\'’ture and density are available only for the 
inner magnetosphere (r s; 15). 

We now present a crude calculation to obtain a value for k. 
The observetl spiralling of the magnetic field lines (Smith et al. 
197h, Northrop et al. 19?4) has been frequently interpreted in terms 
of a radial outflow model (see e.g. Kennel and Coroniti 1975). Tiie 
spiral angle con be related to an outflow velocity very much in the 
same way that the solar wind garden-hose angle is related to the 
solar wind speed. Since the spiralling angle Is proportional to r 


(JorUil’U ) ('i aL., VVf’‘>) t.h« mitl’Inw velocity nhould be conntant. In 
thin can » the AllVen Much mu'bt?r of the flow in 

M. ^ ■' ^1 >. lo* « 0.4 

A y “07^ 

r r 

where we havo used the estimate of Kennel and Coroniti 
V 10 ' km/n. We require that in of order vmlty. (If were 
much luv^er n nhoek woulil fonn and decelerate the flow to 
nub-AliVenie velociticB.) Then k and the density decreonen 

an 


N ' . h X 10' 
o ? — 


At J.4 Hj t!ie dennity would bi» about t cm which is not too different 
from the value Lciu” * quoted by Frank et al. 
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Figure la. The weurped current sheet model (after Smith et al, 
[197^3 and Hill et al. [1974]). 

Figure lb. The rigid current sheet model (after Van Allen et al. 
[1974]), 

Figure 2a. The angle between the model perturbation field 
? = ^2 th(z/D)# 
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and the rotational axis R. The dashed line represents the 
observations of Smith et al. (1974). 

Figure 2b, The angle between the model perturbation field and the 
magnetic dipole axis M. The dashed line represents the ob- 
servations of Smith et al, (1974). 

Figure 3. The magnitude of the model magnetic field strength 
B = along the outuound trajectory of Pioneer 10. 


The dots represent 1 hour averages published by Smith et al. 
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